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ABSTRACT 



Pointing error standard deviations for two theodolites, 
the Wild T-2 and Odom Aztrac, were determined under condi- 
tions closely approximating those of range-azimuth or 
azimuth-azimuth hydrographic surveys. Pointing errors found 
for both instruments were about 1.3 meters, and were inde- 
pendent of distance. No statistical difference between the 
errors of the two instruments was found. The accuracy of 
the interpolation methods used by the National Ocean Service 
(NOS) for range-azimuth positioning were investigated, and 
an average inverse distance of about 2.5 meters was observed 
between interpolated positions and corresponding observed 
positions. The overall range-azimuth position errors of the 
two theodolites were then compared to positioning standards 
of NOS and the International Hydrographic Organization, 
using assumed ranging standard deviations of 1.0 and 3.0 
meters. Both instruments met all standards except the NOS 
range-azimuth standard for 1:5,000 scale surveys. 
Interpolated positions may fail to meet more of the stan- 
dards because of additional inherent error. 



4 



TABLE OP CONTENTS 



I. INTRODUCTION 9 

A. THE RANGE-AZIMUTH POSITIONING METHOD 9 

B. HYDROGRAPHIC POSITION ERROR STANDARDS .... 12 

C. OBJECTIVES 14 

II. ERROR INDICES AND R A NGE-AZI MU T H GEOMETRY 20 

A. DEFINITIONS 20 

1. E landers 22 

2. Systematic Errors 23 

3. Random Errors 24 

B. TWO-DIMENSIONAL ERROR FIGURES 26 

1. Concentric and Eccentric Geometry .... 27 

2. The Error Ellipse 29 

3. Root M«an Squar® Distance 33 

4. Circular Standard Error 35 

C. THE ERROR OF AN INTERPOLATED FIX 37 

1. Interpolation Algorithms 37 

2. Error Propagation 39 

III. EXPERIMENT DESIGN AND IMPLEMENTATION 43 

A. FIELD WORK 43 

3. THEODOLITE POINTING ERROR 47 

C. INTERPOLATION ALGORITHM EVALUATION 49 

D. CHOICE OF EXPERIMENTAL CONDITIONS 51 

IV. RESULTS AND DATA ANALYSIS 56 

A. DATA PROCESSING SYSTEM 56 

3. POINTING ERROR DETERMINATION 56 

C. INTERPOLATION EVALUATION 64 

D. ANALYSIS OF FACTORS AFFECTING THE RESULTS . . 66 



5 



2. APPLICATION TO POSITION ERROR STANDARDS ... 72 



V. CONCLUSIONS AND R ECD MM ENDATIO N S 79 

APPENDIX A: MEANS AND STANDARD DEVIATIONS OF 

ACQUIRED DATA SETS 85 

APPENDIX B: GEODETIC POSITION OF HORIZONTAL CONTROL 

STATIONS 91 

APPENDIX C: EMPIRICAL PROBABILITY DENSITY FUNCTION 

PLOTS 92 

BIBLIOGRAPHY 98 

INITIAL DISTRIBUTION LIST 101 



6 



LIST OF TABLES 



I. 

II. 

III. 

17. 

V. 

VI. 

VII. 
VIII . 

IX. 

X. 

XI. 

XII. 
XIII . 
XIV. 



Circular Error Formulae 

Data Acquisition Sequence of Events 

Pointing Error Standard Deviation (pooled 

estimates) 

Summary of ANCVA Results at 95% Confidence . . . . 

Results of Interpolation Evaluation 

Corrected Aztrac Standard Deviation 

Experiment Precision and Theodolite Error . . . . 

Position Standards Comparison 

3000 Meter Arc 

1500 Meter Arc 

1000 Meter Arc 

700 Meter Arc 

500 Meter Arc 

300 Meter Arc 



36 
4 6 



53 



63 

65 

70 

72 

78 

35 

86 

37 

88 

89 

90 



7 



LIST OP FIGURES 



1.1 Illustration of Range- Azimuth Positioning ... 10 

1.2 Illustration of Aztrac Shore Station 17 

1.3 The Aztrac Transmitting Unit 19 

2.1 The Normal Probability Curve 25 

2.2 Eccentric Range-Azimuth Geometry 27 

2.3 Concentric Range-A zimuth Geometry 28 

2.4 Eccentric Error Compensation 29 

2.5 Error Ellipse and dy-^ s 31 

2.6 A Range-Azimuth Position 32 

2.7 Variation of d Probability 34 

2.8 Example of Angular Interpolation 38 

3.1 Sketch of the Survey Area 44 

3.2 Lines of Position Observed to the Vessel .... 45 

3.3 Uncertainty of an Observed Error 43 

3.4 Interpolation of Angle Only 50 

3.5 Interpolation of Angle and Distance 50 

4.1 One-Dimensional Difference Between Means .... 66 

4.2 Two-Dimensional Difference Between Means .... 67 

4.3 Probability Density vs. Error 68 

4.4 Results Compared to dmi 74 

4.5 Results Compared to 90% Probability 76 

C. 1 Probability Density Plot: 300 m arc 92 

C. 2 Probability Density Plot: 500 m arc 93 

C.3 Probability Density Plot: 700 m arc 94 

C. 4 Probability Density Plot: 1 000 m arc 95 

C. 5 Probability Density Plot: 1500 m arc 96 

C. 6 Probability Density Plot: 3000 m arc 97 



8 



I. INTR ODUC TION 



A. THE RANGE- AZIMUTH POSITIONING METHOD 

The fundamental purpose of a hydrographic survey is 
defined by Ingham (1 574) as being to "depict the relief of 
the seabed, including all features, natural and manmade, and 
to indicate the nature of the seabed in a manner similar to 
the topographic map of land areas." He goes on to describe 
two factors defining a single point on the seabed: 

(i) "The position of the point in the 
horizontal plane in, for example, latitude 
and longitude, grid co-ordinates or anqles 
and distances from known control points. 

(ii) The depth of the point below the sea 
surface, corrected for the vertical 
distance between the point of measurement 
and water level and for the height of the 
tide above the datum or reference level to 
which depths are to be related." 

Thus the hydrographer must answer the two primary gues- 
tions of "hew deep" and "where" for each of the thousands of 
soundings acquired on every survey. Because every area to 
be surveyed has different geophysical characteristics and 
levels of use, the hydrographer must possess a suite of 
tools and techniques to accomplish each survey. A survey of 
a large metropolitan harbor requires different equipment and 
measurement precision than one for a deep ocean area. 

Only the first of Ingham's two factors cited above is 
considered, and it is further narrowed in scope to techni- 
ques used in the most precise surveys. Such a survey might 
be one of a winding, narrow river carrying deep draft 
vessels, or perhaps a very large scale survey of an inner 
harbor. Eoth areas require the highest positioning accuracy 
and a minimum of shore control stations. 
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Any method of positioning employs the intersection of 
lines of position (LOP's) to construct a fix. Although 
advanced methods may use multiple L3P's, traditional hydrog- 
raphy uses the simple intersection of two lines to fix the 
vessel's position. The vessel is located somewhere along 
each of two lines of position, and the only point satisfying 
these conditions is the intersection of the lines. 

The error associated with one of the simplest posi- 
tioning methods, that of the two LDP range-azimuth fix, 
which is illustrated in Figure 1.1, is analyzed. 



Figure 1. 1 Illustration of Range- Azimuth Positioning. 

Also called the rho/theta method, range-azimuth positioning 
consists of the observation of a distance and an azimuth to 
a vessel from either one or two known locations [Ombach, 

1976 ]. An example of this method is the use of radar aboard 
ship. A relative position for a radar contact is determined 
by observing a radar range and azimuth, or a radar range and 
visual azimuth, to a contact. The two lines of position 
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Hy droplet System automatically records a distar.ee measure- 
ment for each sounding. Azimuths ace not telemetered to the 
vessel but are relayed over voice radio and are manually 
entered into the computer system. Since the maximum data 
rate is about two angles per minute, the interpolation of 
angles for sounding positions between fixes is necessary. 

Recently a digital theodolite, the Odom Aztrac, has been 
developed which can record and telemeter angles with great 
speed -- up to ten angles per second [Odom Offshore Surveys, 
Inc., 1982], A computer system aboard the survey vessel can 
thus record and plot an observed position for each sounding. 
This rapid position fixing, combined with a computer's 
ability to provide cross-track error indications to the 
helmsman, enables the hydrographer to systematically cover a 
survey area with maximum efficiency by running straight and 
parallel sounding lines. 

The Aztrac system is still considered a semiautemated 
system because an observer is required to manually track the 
vessel with the theodolite. Two fully automated range- 
azimuth systems which feature fully automatic tracking have 
been devolcped. One is the Polarfix system developed by 
Krupp-Atlas Elektronik in Germany [Smith, 1983], and the 
ether is the Artemis system developed by Christiaan 
Huygenslaborator ium in the Netherlands [Newell, 1981]. 

B. HYDR0GB1PBIC POSITION ERROR STANDARDS 

Historically, most national hydrographic organizations, 
as well as the International Hydrographic Organization 
(IHO) , have used linear plotting error at the scale of the 
survey to be the standard for sounding position accuracy. 
Prior to 1982, the standards recommended by IHO [IHO, 1968] 
were : 



"The indicated repeatability of a fix 
(accuracy of location referred to shore 
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control) in the operating area, whether 
observed by visual or electronic methods, 
combined with plotting error, shall seldom 
exceed 1.5 mm (0.05 in) at the scale of 
the survey." 



The IHO recently published new recommendations for error 
standards [IHO, 1982] which are: 

"... any probable error, measured relative 
to. shore control, shall seldom exceed 
twice the minimum plottable error at the 
scale of the survey (normally 1.0 mm on 
paper) . ” 

Neither of the IHO standards make any reference as to 
what probability level they apply. Munson (1977) inter- 
preted the words ”shall seldom exceed" in the above state- 
ments to mean "less than 10% of the time", which seems 
reasonable. The 1982 IHO standard is somewhat confusing due 
to its use of both the terms "seldom exceed" and "probable 
error". The latter term is associated with a 50% prob- 
ability by most statisticians including Greenwalt (1971). 
However, the author of these standards. Commodore A. H. 
Cooper, SAN (retd) , has stated that he intended no statis- 
tical significance to the term "probable error" [Wallace, 
1983 ]. 

The NOS has not yet incorporated the latest IHO stan- 
dards, but such action is being considered in some form 
[Wallace, 1983]. Current NOS standards have been developed 
to ensure that "accuracies attained for all hydrographic 
surveys conducted by NOS shall equal or exceed the specifi- 
cations" of the 1968 IHO standards [Umbach, 1 976 ]. Unlike 
the international standards, the NOS standards for all elec- 
tronic positioning systems use the concept of root mean 
square error (d rrnS or rmse) , which has a somewhat variable 
probability of between 68.3 and 63.2 percent. The NOS stan- 
dards for fully visual and for hybrid (combination elec- 
tronic and visual) positioning have no explicit reference to 
probability. 
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Specific operational standards for range- azimuth posi- 
tioning have been neglected by many hydrographic organiza- 
tions. However, NOS [Umbach, 1976] requires the following 
observational procedures be followed for all range-azimuth 
positions. 



"Objects sighted on should be at lease 500 
m from the theodolite... the azimuth 
check should not exceed one minute of 
arc... observed azimuths or directions to 
the sounding vessel for a position fix 
shall be read to the nearest 1 min of arc 
or better if necessarv to produce a posi- 
tional accuracy of 0.5 mm at the scale of 
the sur v ey . " 

Since the range-azimuth method is classified as a hybrid 
positioning system, it is not referenced to any particular 
probability, but a reasonable assumption may be made that 
the d riwS concept also applies in this case. 

The U.S. Naval Oceanographic Office (NAVOCEANO) also 
requires that its surveys meet the standards of the NOS 
Hydrographic Manual. The Army Corps of Engineers presently 
have no formal positioning requirements that must be met by 
all districts, although draft specifications are being 
written at this time [Hart, 1983]. The range-azimuth tech- 
nique and its applicability to Corps of Engineers surveys is 
discussed in Hart ( 1 577) . No specific requirements for 
range-azimuth positioning could be found for either the 
Canadian Hydrographic Service or the British Hydrographic 
Service. Palikaris (1983) also reports no published stan- 
dards for these organizations. 

C. OBJECTIVES 

All position error standards using an explicit prob- 
ability are based on the idea that an observation is a 
normally distributed random variable with zero mean and 
standard deviation <r . These standards require a value for 
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confusion when tracking a fast moving vessel. This is 
because the observer sees an image of z he survey boat upside 
down and moving apparently in the opposite direction from 
its actual movement. Another disadvantage is that the mech- 
anism for reading the horizontal circle of a T-2 is more 
complicated than desired for the rapid observations neces- 
sary in hydrographic survey work. The operator is required 
to step tracking the vessel, remove his eye from the tele- 
scope, and use an auxiliary eyepiece to read the angle. 

There is a practical limit of about 30 seconds to the speed 
with which successive angles can be observed and read. The 
particular unit tasted was serial number 30504. 

The Odom Aztrac theodolite is a semi-automated, 
line-of-sight angle measuring system. The Aztrac system 
consists of a Wild T-16 theodolite (serial number 2534880 
was tested) which was modified for infinite tangent drive 
and to provide angular information in a digital format. The 
shore unit decodes the observed angle, determines the direc- 
tion of rotation of the instrument and displays the angle on 
its front panel. The angle is then converted to binary 
coded decimal (3CD) format and used to frequency shift key 
(FSK ) an FM transmitter to link the data with the survey 
vessel. The data is transmitted at the rate of 10 angles 
per second. On board the vessel the Aztrac receiver 
converts the received data to parallel form and displays if 
on the front panel for manual recording. For automated 
recording or processing by onboard computer a serial data 
ouout is provided. The Aztrac has an angular resolution of 
0.01 degree (36 arc seconds) (Odom, 1983 ]. 

With the notable exception of angular resolution, this 
theodolite is more appropriate for range-azimuth hydrography 
than the Wild T-2. It has an erect image and infinite 
tangent screw which allow the vessel to be constantly 
tracked. Its digital output requires no action on the 
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Figure 1.2 Illustration of Aztrac Shore Station, 
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operator's part to record observed angles. An additional 
positive feature of the digital output is the abiliry to 
rapidly zero the instrument on an initial pointing. This is 
done by pressing a manual reset button on the instrument 
control panel. Figures 1.2 and 1.3 show photographs of the 
Aztrac equipment. Figure 1.2 illustrates a typical shore 
station. The observer is adjusting the Aztrac instrument, 
and the transmitter unit is on the ground to the right of 
the Aztrac tripod. The distance measuring equipment is 
mounted on a tripod behind the observer, and the Aztrac 
transmitting antenna is at the top of a pole on the extreme 
right of the photograph. Figure 1.3 shows the Aztrac trans- 
mitter unit, with its digital angular display in hundredths 
of a degree. Both photographs were provided by Odom 
Offshore Surveys, Inc. 

A second objective of this paper is to evaluate the 
interpolation methods used by the MOS for range-azimuth 
work. The availability of the digital theodolite, with its 
direct measurement of all positions, enabled- a comparison of 
observed and interpolated fixes to be made. This made 
possible an estimate of whether interpolated positions meet 
required accuracy standards , and whether direct measurement 
of all positions is needed. The estimate presented here was 
not made statistically rigorous so that the thesis could be 
kept to a manageable size. More theoretical statistical 
work is needed to fully reduce the interpolation data. 

The final objective of this investigation is to compare 
the position errors of these two instruments with the 
various position error standards discussed in section 3 of 
this chapter. The conclusions resulting from this objective 
will assist the hydrographar to select equipment and oper- 
ating conditions that meet required position error 
standards. 
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Figure 1. 3 



The Aztrac Transmitting Unit. 
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II. ERROR INDICES AND RANGE- AZIMUTH GEOMETRY 

The method of r ange- azi mut h positioning is usually 
selected for large scale surveys because of its simplicity 
and accuracy. This is the result of both angle and distance 
measurement devices being co-locatsd, with the intersection 
of the two lines of position always being ninety degrees. 

In practice, however, the co-location of both instruments is 
often not achieved. The result is an eccentric geometry for 
the position fix. This section will analyze the geometry of 
both eccentric and ccncentric fixes. Position error indices 
in common use will be reviewed and analyzed for the special 
cases of range-azimuth methods, and the error of an 
interpolated fix will be derived. 

A. DEFINITIONS 

Although a complete and general treatment of error 
theory will not be presented, some basic definitions are 
necessary to understand the data analysis presented here. 

The ideas in this section are included in many basic statis- 
tics textbooks, and were specifically drawn from Wonnacott 
(1977), Bowditch (1977), Heinzen (1977), Kaplan (1980), and 
Davis (1981) . 

Error may be defined as "the difference between a 
specific value and the correct or standard value" [Bowditch, 
1977], or as "the difference between a given measurement and 
the "true" or "exact" value of the measured quantity" 

[Davis, 1981]. Mathematically it can be defined as: 

e = x. - T (2.1) 



20 



where e is ths error, x,; is an observation, and T is the 
"correct" or "true" value. The word error implies that 
there is a known true value for a quantity, with which a 
measurement may be compared to find the "error" associated 
with that measurement. Since the true value of a measured 
quantity is rarely known, the rerm "error" is not precisely 
correct. Davis (1981) states that it is more appropriate to 
speak of the theory of observations rather than the theory 
of errors, but it can be shown that the difference between 
the two is largely one of semantics. 

A single measurement of a particular quantity may be 
considered sufficient for many purposes, even if it is known 
that additional measurements will probably be slightly 
different than the first. If the quantity to be measured is 
of sufficient importance, then multiple measurements are 
made and the sample mean, "X, is used. Each of these 
multiple measurements can be a considered numerical value 
for a random variable. A random variable is one that takes 
on a range of possible values, each associated with a 
particular probablility. 

The sample mean may be expressed mathematically by equa- 
tion 2.2 [Wcnnacott, 1977]: 



the population mean ^ « . The sample mean is always an esti- 
mate of the population mean, which is never directly 
computed. This leads to the concept of the residual, v, 
which is the difference between ths estimate X of the popu- 
lation mean and the observation x^ . This is shown in 
equation 2.3. 




( 2 . 2 ) 



where n is the sample size. If the sample size were 
increased without limit (n — > o°) , equation 2.2 would give 
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V 



X 



(2.3) 






The residual is computationally the negative of the 
error. Nevertheless, equation 2.3 is more appropriate 
because it uses an estimate, Y, of the unknowable population 
mean . The presence of X in equations 2.3 implies that 
multiple measurements have been made, and allows a partic- 
ular confidence to be assigned to the estimate of 
depending on the number of such measurements. Because the 
word error is still used in much of the hydrographic profes- 
sion, it will be used interchangeably in this paper with the 
term residual. It is important, however, to understand that 
the concept of the residual, whatever its name may be, is 
fundamental to any measurement operation. 

Errors are classically divided into three groups: blun- 

ders, systematic error, and random errors [Greenwalt, 1962]. 
Eowditch (1977) and Davis (1981) do not classify errors as 
including blunders, but like the term error itself, the 
distinction is largely a semantic one. Ideally blunders and 
systematic errors are completely eliminated from the data. 
The most precise measurements reduce random error as much as 
possible, but it can never be completely eliminated. 

1 . Blu nde rs 

Blunders are mistakes that are "usually gross in 
magnitude compared to the other two types of errors" [Davis, 
1981], and are most often caused by carelessness on the part 
of the observer, or by grossly malfunctioning observing 
equipment. They are usually detected and eliminated by 
procedural checks during the data acquistion process. The 
recognition of a blunder is not always easy, since a blunder 
"may have any magnitude, and may be positive or negative" 
[Bowditch, 1977], 
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2 • Systematic Error s 

Systematic errors are defined by Davis (1981) as 
those that occur "according to a system which, if known, 
always be expressed by mathematical formulation." This 
mathematical model results in correctors that are applied to 
all measurements obtained, thus eliminating the systematic 
errors from the observations. The model may be as simple as 
a constant corrector subtracted from lengths obtained with a 
steel tape, or it may be as complicated as modelling the 
effects of atmospheric refraction on electronic distance 
measuring equipment. 

If the systematic error is such that it cannot be 
modelled, it is then estimated by a process known as cali- 
bration. Kaplan (1 980) defines calibration as the process 
of comparing the measuring instrument against a "known" 
standard. The word "known" is usually operationally defined 
as a measurement operation or instrument rhat is much more 
accurate than the one being calibrated. The difference 
between the observed and standard value is used as an esti- 
mate of the total effect of all systematic errors present. 
This process is very close to the classical concept of 
"errors" presented above, and is entirely proper for use in 
the correction of systematic errors [Davis, 1981]. Of 
course, one must be careful to apply the corrector only tc 
those measurments made under the same conditions as the 
calibration. 

A systematic error found in theodolite or sextant 
observations is known as the personal error of the observer 
[Mueller, 1969], [Bowditch, 1977]. This type of error is 
rarely quantified for hydro gaphic applications, but never- 
theless it does exist. The observer must rely on the senses 
of hearing and vision to make measurements, which vary 
between individuals as well as with time in one individual. 
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Some personal errors are constant and seme are erratic 
[Davis, 1981], These errors are minimized by training and 
standardizing observational procedures. The best way to 
eliminate personal error is by the use of completely 
automated observation equipment. 

3. Random E rr ors 



magnitude or sign' 1 , and are "governed by the laws of prob- 
ability" [Bowditch, 1977]. If one assumes that all blunders 
and systematic errors have been removed from the observa- 
tions, the remaining values can be regarded as sample values 
for a random variable. As noted earlier, a random variable 
can take on a range of values, each associated with a 
particular probability. A random error has high probability 
of being close to the population mean , , and a low prob- 
ability of being very much different than [Greenvalt, 

1962 ]. 

A probability density function expresses the rela- 
tion between a value for a random variable and the prob- 
ability of its occurrence. Hydrographic survey measurements 
often use the normal or Gaussian probability density func- 
tion. A concise explanation of this function is given in 
Greenwalt (1962) and Kaplan (1980). The function itself is 
given as equation 2.4, where p (v) is the probability of the 
occurrence of a particular residual v, and <T is the popula- 
tion variance which is approximated by the sample variance, 
S', given by equation 2.5. 



"Random errors are chance errors, unpredictable in 



I V 




(2.4) 



S 







(2.5) 
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To find the probability of a residual falling 
between two residuals v, and v z , equation 2.4 must be inte- 
grated over the interval v, to v z . This corresponds to the 
area under the gaussian curve between “hose two points, as 
is shewn in Figure 2.1. If p(v) were integrated from -1<rto 
+ 1<T, the area under the curve would be 68. 27% of the total 
area. This means that there is a 68.27% probability of a 
particular residual falling between plus or minus one stan- 
dard deviation of the mean, where the standard deviation, <T , 
is defined as the square root of the variance given in 
equation 2.5. 

< 1 



p(y) 




Figure 2.1 The normal Probability Curve. 

The combined effect of blunders, systematic errors 
and random errors can now be seen in overview. If it is 
assumed that blunders and systematic errors have been 
completely eliminated from a set of observations, there 
remain only random errors. If the sample size is large 
enough, then the sample mean and variance are good 
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approximations of the population mean and variance. That 
is, there is a 68. 2755 probability that any future measure- 
ments made under the same conditions will not fall farther 
than plus or minus one standard deviation from the mean. 

B. TKO-DIHENSIONAL EBHOR FIGUHES 

The two-dimensional error of a position, as applied to 
the special case of range-azimuth fixes, must next be exam- 
ined. Two figures, the ellipse and the circle, are used to 
characterize two-dimensional error. The "error diamond" is 
also sometimes used but has no statistical significance 
[Thomson, 1977]. The error ellipse is discussed first since 
it is the most general index of error. Another is root mean 
square distance (d r ^ s ) , also known as root mean square error 
[Bowditch, 1977 ], which is the radius of a circular figure 
commonly used in hydrography. It is the error index used 
for NOS positioning standards [Umbach, 1976]. A second 
circular figure, known as circular standard error, is also 
examined briefly because of the ease of converting it to 
circular figures which have different probabilities. 

For clarification of the issues involved in this 
section, the following assumptions are made. 

(i) Only random errors are considered. 

(ii) Errors associated with each LOP are normally 
distributed. 

(iii) Errors are independent. 

(iv) Errors are limited to the two-dimensional case. 

(v) LO P • s are straight lines at their point 
of i nt e r sec t io n . 
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These are the same assumptions made by Kaplan (1980), 
except that Kaplan allows systematic error to be considered 
in assumption (i). This appears to be an oversight since 
the remainder of his discussion, from which this section 
draws heavily, considers only random error. 

1 . C oncen t ric and Ecce ntri c Geometry 

Before proceeding further into a discussion of error 
figures, it is necessary to examine the two special cases of 
range-azimuth positioning, that of eccentric and concentric 
geometry. Each is illustrated in Figures 2.2 and 2.3. 




In actual practice the geometry used is often eccen- 
tric, but the concentric assumption is made. This is 
because it is usually difficult to co-locate both theodolite 
and ranging equipment. Hence, they are offset one or two 
meters from each other. It should be noted that this 
assumption will introduce a systematic error in all 



27 



1 




Figure 2.3 Concentric Hange- Azimuth Geometry. 

positions, which must be eliminated before an analysis of 
random errors can be made. The systematic error is often 
ignored (by use of the concentric assumption) if the total 
uncompensated error is within the tolerance of the standards 
being used. An algorithm for eliminating this error is 
shown in Figure 2.4 and given by equation 2.6. This algo- 
rithm assumes than the size of d in Figure 2.4 is small 
compared to r. 



c = d cos p 



( 2 . 6 ) 



where: 

c = the corrector to be applied to r 
r = observed range to the vessel 

d = distance between theodolite and ranging device 
p = the angle between the visual LOP and the 
line connecting the two stations 
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Figure 2.4 Eccentric Error Compensation. 



2. The Error El lips e 



Detailed discussions of the development of the error 
ellipse can be found in many references, especially in 
Greenwalt (1962) and in Burt (1966). This paper will only 
present enough background to apply the error ellipse concept 
to two LOP range-azimuth positioning. The error ellipse 
formed when multiple LOP observations are made is not 
considered here. 

A range-azimuth position is formed by the intersec- 
tion of two LOP's, each having an associated standard devia- 
tion. By applying the two-dimensional normal distribution 
to the errors, elliptical contours of egual probability 
density are formed. The contours cenrer on the intersection 
point of the lines of position. This is illustrated in 
Figure 2.5, and shown mathematically by 



P (v a ,v b ) 




(2.7) 
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where : 

\ = residual (error) in the direction of 
\ - residual (error) in the direction of <J£ 
<3* = length of the semi-major axis 
<T b = length of the semi-minor axis 

It can be shown that 




( 2 . 8 ) 



wh e r e 



K 



- Z JU 



p(v«, v b) <r a <r b 2ir 



( 2 . 9 ) 



'•For values of p (v , v ) from 0 to oo , a family of equal 
probability density ellipses are formed with axes K and 
K(T t " [Greenwalt, 196 2], The probability density function i 
equation 2.7, when integrated over a particular area, 
becomes the probability distribution function. This yields 
the probability that the residuals v^ and v to will occur 
simultaneously within that region. This probability distri 
bution function of an ellipse is given by 

-K/z 

P(v A ,v b ) = 1 - e (2. 10) 



The solution of equation 2.8 for different values of K 
yields different probabilities. For example, for 39.35% 
probability, the axes of the ellipse are 1.00 0 ck and 1.000 
C j" b [ Greenwalt, 1 962 ]. In other words, a one-sigma error 
ellipse around a measured position indicates a 39.35% prob- 
ability that the position is actually within that ellipse. 

It is seen in Figure 2.5 “hat the standard devia- 
tions^ and <3[ of the measured LOP’s are not the same as the 
standard deviations and c£ of the error ellipse. A 
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Figure 2.5 Error Ellipse and d rmS 



coordinate transformation is required to obtain them. This 
transformation is found in Heinzer. ( 1 977) and will not be 
discussed here. Results of the transformation are presented 
in equations 2.11 and 2.12, from Bouditch (1977). 



<T* = 



< = 



2 sirf/3 



2 sin/S L 



<J t + <3; + ~\J cr‘)-4siryS(r>^ 

0 ; - <r‘ - ^(<^<Q-4sin/3<r>* 



( 2 . 11 ) 



( 2 . 12 ) 
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where : 



<%x = length of the semi-major axis 
Gb = length of the semi-minor axis 
°i = standard deviation of the range LOP 
Ok = standard deviation of the angle LOP 
when converted to distance units 

( 2 = angle of intersection of LOP’s 

For range-azimuth positioning, 07 and are tot 
equal. (T, is the error in distance measurement, and it is 
dependent on the equipment used for ranging. A diagram of a 
range-azimuth position is shown in Figure 2.6. 




Figure 2.6 A Range-Azimuth Position. 
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